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Abstract: 

We propose a simple method to calculate the pion form factor at not very large momentum 
transfers, which combines the technique of the QCD sum rules with the description of the 
pion in terms of the set of wave functions of increasing twist. This approach allows one to 
calculate the soft (end point) contribution to the form factor in a largely model-independent 
way. Our results confirm existing expectations that the soft contribution remains important 
at least up to the momentum transfers of order 10 GeV 2 , and suggest that it comes from 
the region of relatively small transverse separations of order 1 GeV -1 . 
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1. Hard exclusive processes in QCD |jj are attracting continuous interest for already 
two decades. In difference to inclusive reactions, like the deep inelastic scattering, exclu- 
sive processes are selective to the partonic content of the participating hadron. Apriory, 
one could consider two different possibilities to transfer a large momentum to a hadron. 
One is the so-called Feynman mechanism, in which large momentum transfer selects the 
configuration, in which one parton carries almost all the momentum of the hadron. The 
transverse size of the remaining "soft" cloud remains in this case arbitrary. The second is 
the hard rescattering mechanism, in which large momentum transfer selects configurations 
with a small transverse size and a minimal number of Fock constituents. In this case the 
momentum fraction carried by the interacting parton (quark) remains an average one. The 
hard rescattering mechanism involves a hard gluon exchange, and can be written in the 
factorized form [Q. 

It has been proved || that the hard rescattering mechanism is the leading one at asymp- 
totically large Q 2 , and yields the pion form factor 

2 

(1) 

Here 4> n {u) is the pion wave function of the leading twist, which describes the distribution 
of the valence pion constituents in the longitudinal momentum. Note that convergence of 
the integral in (|l|) requires that the pion wave function decreases at u — > 1, and the crucial 
point in establishing the formula in ([I]) was the proof [§] that at asymptotically large Q 2 
the wave function is given by the simple formula 

4">(u) = 6u(l-u). (2) 

The results in ([{]), (0) belong to the most important and most rigorously proved statements 
in QCD. 

At large but finite momentum transfers there might be a number of corrections to the 
hard rescattering formula in (P, and till now there was a much more moderate progress in 
understanding whether available energies could be treated as asymptotic ones. An attempt 
to describe the data for the pion form factor starting from Q 2 > 3 GeV 2 by the contribution 
of hard rescattering alone, implies that the low energy pion wave function must be very 
different from its asymptotic form, an issue which has been put to fore and studied in detail 
by Chernyak and Zhitnitsky ||. Using the QCD sum rule approach they have shown that 
the pion wave function at a low scale is wider than the asymptotic one, and proposed a 
model 

^' z \u,fx ~ 500 Me V) = 30u(l - u)(2u - l) 2 , (3) 

which has a peculiar "humped" profile, with a zero in the middle point, corresponding 
to the symmetric configuration where the quark and the antiquark carry equal momenta. 
However, there exists a number of arguments that force to doubt the assertion about the 
dominance of the hard-scattering mechanism in the region of available Q 2 ~ 1 — 10 GeV 2 . 
The well known point of Isgur and Llewellyn-Smith 0] is that the wave functions of the type 
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suggested in || strongly emphasize the contribution of the end-point region of large u in (|IJ), 
where the virtuality of the gluon is not enough to justify the perturbative treatment. This 
contribution of large u, 1 — u ~ m 2 /Q 2 (where m 2 is a certain hadronic scale) corresponds to 
the Feynman mechanism to transfer the large momentum, and should be treated separately. 
At present, there is an increasing evidence that this contribution is numerically important 
up to very high Q 2 , although it is down by an extra power 1/Q 2 in the asymptotics. 

Using the QCD sum rule approach |j] it has been shown || [7J that the pion form factor at 
Q 2 ~ 1 — 2 GeV 2 is practically saturated by the Feynman- type contribution. Unfortunately, 
this method in its standard form cannot be applied for higher values of Q 2 , since it involves 
the expansion in the contributions of vacuum condensates, which coefficients appear to be 
enhanced by increasing powers of Q 2 . Thus, at sufficiently large momentum transfers the 
expansion breaks down. In attempt to cure this problem Radyushkin and collaborators 
H |9| have suggested to resum the series of power corrections by introducing the nonlocal 
extension of the concept of the vacuum condensates, which takes into account the final 
correlation length in the QCD vacuum. Results of || indicate the dominating role of the 
soft contribution at least up to Q 2 ~ 10 GeV 2 . 

The approach of || clearly demonstrates the origin of difficulties in the standard QCD 
sum rule calculations, but may receive objections concerning its theoretical accuracy, since 
not all high-order contributions can consistently be taken into account in this way, and 
also the parametrization of the nonlocal condensates is essentially model-dependent. In 
this letter we suggest an alternative approach to the calculation of the pion form factor 
at not very large values of Q 2 , which seems to be free from both the difficulties of the 
standard QCD sum rule method, and the ambiguities involved in its extension in || [|. 
Our method essentially follows earlier works flQfl-|15||, where the QCD sum rule approach 
has been modified to incorporate the operator product expansion in powers of the deviation 
from the light-cone (in contrast to the short distance expansion in ||). 

2. The idea is to combine the standard technique for the study of hard exclusive 
processes and the QCD sum rule method. To this purpose, we consider the correlation 
function 

T» u (p,q) = i J dxexp(tqx)(0\T{3l(0)f™(x)}\n + (p)) , (4) 

where ft = d'j^u and H m = e u u^ u u + edd^d is the electromagnetic current. At large 
Euclidian momenta (p — q) 2 and q 2 this correlation function can be calculated in QCD, in 
the precise analogy with the calculation of the 717*7* form factor (for different virtualities 
of photons), 0. The leading contribution is written in terms of the pion wave function of 
the leading twist 

T^{p,q) = ZifnPuPv I du n ^ L ~l ( 5 ) 

Jo (1 — u)Q z — u(q — p) z 

where is the pion decay constant, Q 2 = —q 2 , and the ellipses stand for the contributions 
of other Lorentz structures. On the other hand, the dispersion relation over (p — q) 2 relates 
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the correlation function (0) to the pion form factor: 



= if^p - q)^ m2 _ _ ^ 2 2F n {q 2 )p v + ... , (6) 

where the dots stand for higher resonances and the continuum contributions. Matching 
between the representations in @ and @ at no-so-large Euclidian —(q — p) 2 ~ 1 GeV 2 , 
we obtain a sum rule for the pion form factor. 

To this end, we use the standard concept of duality, which tells that the pion occupies 
the "region of duality" in the invariant mass of the qq pair, up to a certain threshold 
So ~ 0.7 — 0.8 GeV 2 . Note that the formula in (|5|) can be rewritten as the dispersion 
relation in {p — q) 2 , with s = (1 — u)Q 2 /u being the mass of the intermediate state. To pick 
up the contribution of the pion, we cut the dispersion integral at s = s , which translates 
to a low bound for the integral over u: u min = Q 2 / (s + Q 2 )- Following ||, we also use the 
Borel transformation to convert the power suppression of higher mass contributions in the 
dispersion integral to the exponential suppression 



u 



(1 — u)Q 2 — u(q — p) 2 
1 



m 2 — (q — p)' 




exp 



where M 2 is the new variable (the Borel parameter) . In what follows we put the pion mass 
to zero. 

Thus, the sum rule arises 



F n (Q 2 ) = jf duUu) exp - K j£* 8«- 



s + Q 2 



which should be satisfied at the values of the Borel parameter M 2 of order 1 GeV 2 . The 
pion wave function in (g) should be taken at a low normalization scale, of order of the Borel 
parameter. 

In what follows we shall complement the sum rules in @ by contributions of higher 
twist. Before doing this, and before going over to the numerical analysis, let us study the 
behavior of the sum rule in the limit of large momentum transfers Q 2 — > oo. 

Because of the 0-function, the integration region in (§) is restricted to values 1 — 
u < so/(so + Q 2 ) — > 0. Thus, the form factor is sensitive to the wave function in the 
highly asymmetrical configuration, where the scattered quark carries almost all the pion 
momentum. According to the general analysis in 0, |3| the behavior of the pion wave 
function in this region coincides with the asymptotic behavior in (Q), ^(u) "~ 1 1 — u. 
Thus, asymptotically, the sum rule in yields 

FjeVff^e-" (9) 
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where (p' n (0) = (d/du)<p n (u)\ u=0 . (We have made use of the symmetry (f> n (u) = (j) n (l — u)). 

For the asymptotic wave function (0) (j)' n (0) = 6. For the model suggested by Chernyak 
and Zhitnitsky (|) </C (cz) (0) = 30, and thus it gives rise to a much bigger soft contribution 
to the form factor. This is in agreement with the observation by Isgur and Llewellyn-Smith 
|3J] that "realistic" model wave functions, which allow for the description of the data on pion 
form factor at virtualities of order a few GeV by the hard rescattering mechanism alone, in 
fact greatly enhance the soft contribution coming from the end point region. 

3. Now we are in a position to take into account the corrections to the sum rule 
in (d), induced by the pion wave functions of the nonleading twist 4. These corrections 
correspond to the contributions to the correlation function in (^), which are suppressed by 
an extra power of the large momentum (p — q) 2 or q 2 . Physically, these corrections have two 
sources. One of them is the transverse momentum of quarks in the leading order diagram 
in Fig. 2a. The second source are contributions of higher Fock states in the pion wave 
function — containing gluon fields in addition to the quark and antiquark ones. In fact, 
these two effects are indistinguishable thanks to the equations of motion, which allow one 
to eliminate all transverse degrees of freedom at the cost of introducing the higher Fock 
components. In this letter we cannot give all the details of the calculation, which becomes 
rather technical. A suitable formalism for the operator product expansion of the correlation 
functions beyond the leading twist has been developed in [ |l~6|j . 

The complete contribution of the diagram in Fig. la to the twist-4 accuracy reads 

n (2o) = Itp^U J q du^ - j—^ - An 5 - up)4 + 2u j + ■■■ ,(10) 

where we have introduced the pion wave functions of twist 2 and 4 defined by the matrix 
element [|Tjj 

(0|d(0)7 M 75«(a;)|7r(p)) = if^ f due^^u) + x 2 9l (u) + 0(x A )) 

2 -i 

+ - I due- mpx g 2 {u) + ..., (11) 

px Jo 

and all logarithmic dependence on x 2 is implicitly included in the wave functions. The 
function G 2 in ( |i~0"D is defined as g 2 (u) = —(d/du)G 2 (u). The diagram shown in Fig. lb 
produces an additional contribution (03 — 1 — a% — a 2 ) 



11(26) = 2i/„ 



1 udu 
(q — up) 4 Jo 




where we have followed in the definition of three-particle wave functions of twist 4: 



(0\d(-x) 1 , l5 gG aP (vx)u(x)\n(p))=p^p a xp-ppx a )^uJ Da^a^e-^-^^ 
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+ 



M<W — — ) -p4<W 



f n / Da$±(ati)e 



-ipx(ai—a2+va3) 



.(13) 



px px 

The wave functions ^\\ and are defined similar to (0), with the substitution (jsgGap) 

A systematic study of the higher twist wave functions has been done in the work |Ti 
and makes use of the expansion in representations of the collinear conformal group SO(2,l), 
which is a subgroup of full conformal group acting on the light-cone. The asymptotic wave 
functions are defined as contributions of operators with the lowest conformal spin. The 
set of wave functions suggested in |L7j includes also the corrections corresponding to the 
operators with the next-to-leading conformal spin, which numerical values are calculated by 
the QCD sum rule method. The results for twist 4 wave functions are (hereafter u = 1 — u): 

120e<5 2 (ai — a 2 )aia 2 a 3 , 



*||(o<) 
*±(ai) 

92(u) 
G 2 (u) 



-120<5 2 aia 2 a3 - + e{l — 3a 3 ) 



305 2 (a 1 — a 2 )a^ 
30S 2 a 2 3 (l- a 3 ) 



2e (1 - 2a 3 



+ 2e(l - 2a 3 



(14) 



95 r 

— S 2 u 2 u 2 + eS 2 \uu{2 + IZuu) 
6 L 

+ 10u 3 (2 - 3u + -u 2 ) In w + 10w 3 (2 - 3u + -u 2 ) In u 
5 5 
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S 2 uu(u — u) , 



5 2 ~ 0.2 GeV 2 , e ~ 0.5. 
Adding the higher-twist contributions in ( |T0"D and (|T2| ) to the sum rule in 



F n (Q 2 



+ 



du exp 



uQ 2 



uM 2 Jo 



d(X\ I da 2 



uM 2 



4 2 

M u ) ~ + G 2 (u)) + Jp92{u) 



1 — 2u + «i — Ct 2 / ^ 
+ 5^ -($|| + 2$_ 



>eU- 



we arrive at 
(15) 

Q 2 \ 



s + Q 2 



a 3 a 3 

which presents our final result. Note that the higher-twist contributions are suppressed by 
a power of the Borel parameter M 2 , as expected. 

With the particular expressions ( [15] ) the integrals of the three-particle wave functions 
can be taken analytically, yielding 

/ da\ [ da 2 — — — = — 5 2 uu{l — 2uu) 

Jo Jo " 



a 3 

, 1 - 2u + ax — a 2 . n , . 
dax I da 2 ^ ($|| + 2$j_) 



10 _ 9 / 15 \ 
-2gx{u)-—5 2 uu(l-—uu). (16) 



Note that the term ~ uu in (|T6| ) cancel in the sum rule (ITq), and thus contributions of three- 
particle wave functions turn out to be of order 1/Q 6 , i.e. are suppressed by an additional 
power of 1/Q 2 - However, the twist 4 contributions still survive in the high-Q 2 limit due to 
the contribution of 5% which yields the same asymptotic behavior ~ 1/Q 4 as the leading 
twist contribution in (|S|). 

4. We turn now to the numerical analysis. Apart from the wave functions, the sum 
rule in ( |TED depends on the value of the continuum threshold Sq, and on the Borel parameter 
M 2 . We take s = 0.7 - 0.8 GeV 2 and vary M 2 in the interval 1-2 GeV 2 , which is 
the expected stability region. The results are shown in Fig.2. In Fig. 2a we plot the value 
of Q 2 F n (Q 2 ) as a function of Q 2 for s = 0.7 GeV 2 and s = 0.8 GeV 2 and for different 
choices of the leading twist pion wave function 7r ( - u): asymptotic wave function (Q) and 
the Chernyak-Zhitnitsky model (^). Since this model in fact refers to a substantially lower 
normalization point than the typical value of the Borel parameter in the sum rule, we give 
the results also for the Chernyak-Zhitnitsky wave function rescaled to /x 2 ~ 1 — 2 GeV 2 

<j>£ Z \u^ ~ lGeV) = 6u(l - u)[l + 0.44C 2 3/2 (2m - 1)] . (17) 

The wave function in (|17D corresponds to the value of the second moment ((2m — l) 2 ) = 0.35 
T2j, which is to be compared to ((2m — l) 2 ) = 0.43 for (|3|). We remind that for the 



asymptotical wave function ((2m — l) 2 ) = 0.2. The contribution of wave functions of twist 
4 does not exceed 20%, and these wave functions are not far from their asymptotic form. 



Thus possible inaccuracy in the model wave functions in ({15) does not have any noticeable 
effect. The stability of the sum rule (|T6|) to the choice of the Borel parameter is illustrated 
in Fig. 2b for several values of Q 2 . 

It is seen that the soft contribution to the pion form factor clearly dominates at Q 2 ~ 
1 — 3 GeV 2 and constitutes about 15-30% of the experimental value at Q 2 ~ 10 GeV 2 (for 
the asymptotic wave function). For the Chernyak-Zhitnitsky model, the soft contribution 
increases substantially. 

Within our approach, the hard gluon exchange contribution originates from the radiative 
correction to the contribution of the leading twist, see diagram in Fig. lc. This contribution 
is not restricted to the end-point region, and thus has no power-like 1/Q 2 suppression. Its 
explicit calculation goes beyond the tasks of this letter. As a rough estimate, one can 
use the expression in (|1|), yielding Q 2 F^Q 2 )^^ ~ 0.15 and Q 2 F v (Q 2 ) hSLld ~ 0.3 for the 
asymptotic wave function and the Chernyak-Zhitnitsky model, respectively. One sees that 
the full answer for the pion form factor, given by the sum of the soft and hard contributions, 
is likely to overshoot the data, if one uses the Chernyak-Zhitnitsky model. 

Main lesson to be learnt from our calculation is that the soft contribution to the pion 
form factor decreases very slowly with Q 2 and is important in the whole region of momentum 
transfers, which are available at present. This conclusion is in full agreement with the results 
of H, |9], although our argumentation is different. 

5. The method described above is quite general, and can be applied to different 
form factors as well. As a one more example, we calculate here the soft contribution to the 
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transition form factor 7717). For the |A| = 1 p - meson (p_|_ hereafter) the transition form 
factor is defined as 

(p±(pi)|#>(P»)> = wM^pW 2 ) , (18) 
where Q 2 = — (pi — P2) 2 and e p ± stands for the polarization vector of the p. This process is 
clearly due to non-leading twist effects || as long as it is related with the helicity flipping. 
As a result, the hard rescattering diagram yield the asymptotic behavior F vp (Q 2 ) ~ l/Q 4 
||. As we will show, the soft contribution exhibits the same asymptotic dependence 1/Q 4 , 
and is of the same order, therefore, as the hard contribution (up to the Sudakov suppression, 
which is unlikely to be important at moderate values of Q 2 ). 
We consider the correlation function 



A(p,q) = i J dxe^(0\T{d(0)a^u(m m ( x )}\<P)) 

= (e^xaq^ - £^\aqu)q\Pa'n(q 2 ,(p - qf) + ••• (19) 

which contains a contribution of interest of the p^-meson 



n '^-«» = g ro ;-T-V (20) 

Here /J is the p^-meson decay constant: 

(0|d(0)^ 76 «(x)|7r(p)) = i(efp u - e^)fj . (21) 
Calculation of the diagram in Fig. la yields in this case 

I% 2 , (p - q) 2 ) = -(e u + e d ) Um } f du/^- , (22) 

3(m u + m d )Jo {q-upp 

where <f a ( u ) is the pion wave function of twist 3 JHj 

— i f Tn 2 

(0\d(0)^ vl5 u(x)\7r(p)) = - U * M x v -p v xt:) / due-^ip^u) . (23) 

6(m u + m d ) Jo 



It has been shown in [T7| that the wave function (p a (u) is close to its asymptotic form 
<f a (u) = 6u(l — u). 

It is easy to check that the diagram of Fig. lb does not contribute to the Lorentz structure 
of interest. Then, to the twist 3 accuracy, we arrive at the very simple sum rule 

F„ P (Q) = (e u + e d ) 3U r M2 J q due &—e( u -—^), (24) 

in which we have replaced the factor appearing in the normalization of the wave function 
<f) a (p3|) by the quark condensate (qq) ~ —(250 MeV) 3 . Following 0, |3| we use the values 
So = 1.5 GeV 2 and fT ~ 200 MeV for the continuum threshold in the p-meson channel, and 



the p-meson coupling, respectively. The results are shown in Fig. 3. Using them we obtain 
an estimate for the \I/ — > 7 — > 7r u> decay rate (F nuJ (Q 2 ) = 3F np (Q 2 ) due to the isospin 
symmetry) 

Br _> 7 _> ^ /m _> e + e -) = A( M *F™(M|)) 2 ~ (4 ± 2) x 10~ 4 . (25) 

A relatively large error is due to a poor stability of the sum rule in this case. The number 
in ( f25|) appears to be in good agreement to the experimental number (4.2 ±0.6) x 10~ 4 flI8| 



The contribution to this form factor of the hard rescattering has been calculated by 
Chernyak and Zhitnitsky ||, using the leading twist pion wave function in ([3|), and three- 
particle wave functions of the p-meson of nonleading twist. It has the same functional 
dependence ~ 1/Q 4 , and approximately the same numerical value as the soft contribution 
which we have calculated here. A simple patching them together would yield the branching 
ratio ^ — > tt°uj several times above the data. To our opinion, the hard contribution to this 
decay given in || is strongly overestimated. 

6. In this letter we have suggested a simple method to calculate the pion form 
factor in the region of intermediate momentum transfers, which is essentially a hybrid of 
the standard QCD sum rule approach and the conventional expansion in terms of the pion 
wave functions. Its value is in the possibility to estimate the soft (end point) contribution 
to the form factor in a model independent way, which is a problem of acute interest. The 
main advantage compared to the standard QCD sum rule calculation H 0] is that the 
"light-cone sum rules" suggested in this paper remain well-defined in the limit Q 2 — > 00, 
and is related to the fact that the parameter of the expansion in our sum rules is the 
twist of relevant operators, but not their dimension as in the standard sum rules. In this 
way contributions of various local operators are resummed in the set of wave functions of 
increasing twist, the end-point behavior of which is known from general arguments. In effect, 
explicit factors ~ Q 2 which may appear in the calculation of higher twist contributions will 
be compensated by factors ~ l/Q 2 originating from a more fast decrease of higher twist 
wave functions at u — > 1 compared to the leading twist ones. The physical reason for 
disappearance of the contributions enhanced by powers of Q 2 is in our approach the same 
as in the calculation involving the nonlocal condensates in ||. However, our method is 
practically model- independent. 

Our main result is the calculation of the soft contribution to the form factor, which 
turns out to be large at least up to Q 2 ~ 10 GeV 2 . In agreement to f§] we find that this 
contribution depends strongly on the shape of the pion wave function. Patching together 
the contribution of hard rescattering and the soft contribution, we find that the model by 
Chernyak and Zhitnitsky is likely to overshoot the data. Combining this result with the 



criticism in || 112] , we conclude that there is increasing evidence, coming from different 
calculations, that the true low energy pion wave function is not that much different from 
its asymptotic form, as proposed in J|. 

On the evidence of an impressive calculation of Sudakov-type double-logarithmic cor- 



rections to the contribution of the hard rescattering, it has been argued in [19| that the 
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end-point contribution to the pion form factor is strongly suppressed already at moderate 
Q 2 ~ 10 GeV 2 . Radiative corrections to the correlation function in (f|) can involve loga- 
rithms of the type ln(g 2 /(g — p) 2 ), but not \ii(Q 2 /Aq CD ), since it is IR-protected. Hence 
the corrections to the sum rule can only be accompanied by logs like ln(Q 2 /so), ln(Q 2 /M 2 ) 
which never become large (at moderate momentum transfers). Thus, the Sudakov expo- 
nential suppression is not likely to occur in our sum rules. To our opinion, the significance 
of Sudakov corrections is overestimated in [PS]. The reason is that the effective transverse 



momentum, generated by the Sudakov suppression, should be compared not to the Comp- 
ton wave length of the pion, of order 1/200 MeV, but to the average transverse separation 
of the quark and antiquark in the particular configuration which dominates the soft con- 
tribution to the form factor. Note in this respect, that the quark- antiquark separation in 
a "free" pion does not have any physical meaning beyond the leading twist accuracy, since 
effects of transverse degrees of freedom can be rewritten in terms of higher Fock components 
in the wave function thanks to the equations of motion. For the particular hard process, 
however, the question of relevant transverse distances is well-defined. In the case of the 
correlation function in ([|), the characteristic transverse separation between the quark and 
the antiquark is given by the deviation from the light-cone x 2 ~ (q — up) 2 , as can easily 
be checked by an explicit calculation in light-cone coordinates in the position space. After 
the Borel transformation, (q — up) 2 goes into uM 2 , so that the characteristic transverse 
separations yielding the form factor in Fig. 2 are of the order 

(26) 

Note that l/(uM 2 ) is exactly the expansion parameter in our calculation, which controls 
the size of higher-twist corrections, and in the working region of the sum rule is of order 
1/sq ~ (0.2 — 0.3 fm) 2 . To our opinion, it is this scale rather than Aqcd which should serve 
as the IR cutoff in the calculation in [[nj. Since the average value of u under the integral 
in the sum rule increases slightly with Q 2 , one may speculate that the relevant transverse 
size and the importance of higher twist contributions are slightly decreasing. 

Thus, the soft contribution to the pion form factor comes from configurations with a 
much smaller transverse size than the poin electromagnetic radius ~ 0.65 fm, which is 
dominated by contributions of multiparton states. At distances ~ 0.2 — 0.3 fm the strong 
coupling is already not large, and the application of perturbation theory to the calculation 
of the contribution of the hard gluon exchange can be justified. However, this contribution 
must be complemented by the contribution of Feynman type, coming from the end-point 



region. Our conclusions essentially support the picture described in p0| 
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Captions 



Fig. 1 Leading contributions to the expansion of the correlation function in (f|) in powers 
of the deviation from the light-cone. 

Fig. 2 The soft (end point) contribution to the pion electromagnetic form factor for the 
asymptotic pion wave function and for the Chernyak-Zhitnitsky model as a function 
of Q 2 (a) and in dependence on the Borel parameter M 2 in the sum rule (b). The 
solid and dashed curves in Fig. 2a correspond to the calculation with sq = 0.8 and 
So = 0.7 GeV 2 , respectively, and the value of the Borel parameter M 2 = 1.5 GeV 2 . 
The curves in Fig. 2b are calculated using sq = 0.8 GeV 2 and the asymptotical wave 
function. Among the pairs of curves marked "CZ" the upper ones correspond to the 
calculation using the Chernyak-Zhitnitsky wave function at the scale fi = 500 MeV, 
and the lower ones at [A — 1 GeV, see and (|T7|), respectively. 

Fig. 3 The transition form factor 7p7r as a function of the momentum transfer. 
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Figure 1: Leading contributions to the expansion of the correlation function in (|J) in 
powers of the deviation from the light-cone. 
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Figure 2: The soft (end point) contribution to the pion electromagnetic form factor for 
the asymptotic pion wave function and for the Chernyak-Zhitnitsky model as a function 
of Q 2 (a) and in dependence on the Borel parameter M 2 in the sum rule (b). The solid 
and dashed curves in Fig. 2a correspond to the calculation with s = 0.8 and s = 0.7 
GeV 2 , respectively, and the value of the Borel parameter M 2 = 1.5 GeV 2 . The curves in 
Fig. 2b are calculated using s = 0.8 GeV 2 and the asymptotical wave function. Among 
the pairs of curves marked "CZ" the upper ones correspond to the calculation using the 
Chernyak-Zhitnitsky wave function at the sca3e \x = 500 MeV, and the lower ones at jj, = 1 
GeV, see @ and ([17|) , respectively. 



Figure 3: The transition form factor 'ypir as a function of the momentum transfer. 
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